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^ \ ABSTRACT 

We use Herschel-PACS far-infrared data, combined with previous multi-band information 
and mid-IR spectra, to properly account for the presence of an active nucleus and constrain 

■ its energetic contribution in luminous infrared (IR) sources at z ~ 2. The sample is composed 
I of 24 sources in the GOODS-South field, with typical IR luminosity of 10'^ L©. Data from 

the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN 
emission. 

We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy 

■ and an AGN component. A smooth-torus model for circum-nuclear dust is used to account 

for the direct and re-processed contribution from the AGN. 

We confirm that galaxies with typical L8_iooo^m~lO^^L0 at z~2 are powered predominantly by 
star-formation. An AGN component is present in nine objects (~35% of the sample) at the 3cr 
confidence level, but its contribution to the 8-1000 fim emission accounts for only ~5% of the 
. energy budget. The AGN contribution rises to ~23% over the 5-30 fim range (in agreement 

rS ■ with Spitzer IRS results) and to ~60% over the narrow 2-6 jt/m range. The presence of an AGN 

I is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating 

" " " the presence of significant absorption, i.e. N//~10^^ - 10^** cm"^. An additional source shows 

indications of obscured AGN emission from X-ray data. The comparison between the mid- 
IR-derived X-ray luminosities and those obtained from X-ray data suggests that obscuration 
is likely present also in the remaining six sources that harbour an AGN according to the SED- 
fitting analysis. 
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1 INTRODUCTION 

Ultra-lu minous infrared ga laxies (ULIRGs, defined as having Ljr > 
10'^L^: ISoifer et al.lll987l) are among the most luminous objects 
of the Universe, radiating most of their energy in the infrared 
(IR) band. In the local Universe, ULIRGs are rare objects and, 
although very luminous, acc ount for only ~5% of the total inte- 
grated IR luminosity density dSoifer & Neugebauej|l99lh . Moving 
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to higher redshift, the contribution of ULIRGs to the energy bud- 
get incr eases, as clearly illustrated by infrared lum i nosity function 
studie s jLe Floc'h etallllOOSi : Icaputi et alJboOTh . iMagnelU et alj 
H), deriving the infrared luminosity function up to z~2.3 from 
deep 24 and 70 yum Spitzer data, obtained a contribution of ~17% 
from ULIRGs to the IR luminosity density at z~2. This result has 
been recently confirmed also by the H erschel Science Demon stra- 
tion Phase (SDP) preliminary results; iGruppioni et al.l ( I2OI0I) . by 
deriving the total IR luminosity function up to z ~ 3, estimated that 
ULIRGs account for ~30% to the IR luminosity density at z~2. 
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At high redshift (z~2), a key question regards the nature of the 
sources that power these ultra-luminous objects (i.e., star-formation 
vs. accretion). Locally, thanks to Infrared Space observatory (ISO) 
mid-IR spectroscopy, ULIRGs were proven to be powered predom- 
inantly by star formation in the mid-IR, with the fraction of AGN- 
powered objects increasing with luminosity ( from ~15% at Lm < 
2x10'^ Lo to about 50% at higher luminositv: iLutz et aljfiggsl) . A 
limited percentage (15-20%) of mid-IR light dominated by accre- 
tion processes has also been found using L-ba nd (3-4 ^m) observa - 
tions of local ULIRGs with ~ IO'^Lq bv lRisaliti et"d] bOld) . 
Interestingly, these authors show that, although their sources are 
powered mostly by starburst processes, at least 60% of them con- 
tain an active nucleus. Using accurate optical spectral line diag- 
nostics applied to a samp l e of 7 jim selected luminous sources 
at z -l. lSvmeonidis et alj feoid) found that only 20-30% of their 
objects may host an active nucleus. However, such AGN are never 
bolometrically dominant. In the same paper, a discussion on how 
the low AGN incidence can be partially due to selection effects 
(i.e. the 70 yum band sampling the starburst 50 yum rest-frame far-IR 
bump) is also presented. 

At z ~ 2, our understanding of the AGN content in ULIRGs 
is much more uncertain. On the one hand, sources selected with 
a 24 yum flux density i 0.9-1 mjy mainly sample the bright end 
of the ULIRG population. The analysis of the IRS spectra of 
a color-selected sample of sources with 824^,,, ~ 0.9 mJy and at 
z > 1.5 (L/,;~7x10'^Lq) sh ows that the major ity of these sources 
are AGN-dominated (~75%: ISai"ina et alj2008h . On the other hand, 
sub-millimeter selected galaxies (SMGs) at z~2, falling in the 
ULIRG regime, appear largely starb urst-dominated objects (e.g., 
IValiante et al.ll20oi|Pope et alj2008l) . 

In this paper, we aim at providing a better understanding of 
luminous sources at z ~2 with typical IR luminosities of IO'^Lq, 
i.e. sampling the knee of the IR luminosity function instead of 
the br ight end. We will re-analyse the sample from Fadda et al. 
( I2OIO . hereafter FIO), where an estimate of the AGN contribution 
in ULIRGs has already been presented, based on Spitzer IRS data, 
Chandra 2 Ms observations, optical/near-IR multi-band properties, 
and ACS morphological properties. In this work, we will extend the 
analysis to far-IR data, obtained recently by the Herschel satellite 
as part of the gu aranteed survey "PACS Evolutionary Probe" (PEP; 
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201 ih . and to the recently published Chandra 4 Ms data 



201 ih . 



We note that an analysis of the FlO sam ple (including the new 
far-IR data) has been already presented by ' Nordon et al.l (1201 2|) 
[hereafter N12], where the general properties of the mid-to-far IR 
spectral energy distributions (SEDs) of 0.7< z <2.5 galaxies have 
been investigated. The authors confirm the early Herschel results 
(e.g.. lElbaz et"al]|2010l) , i.e. the star-formation rates at z ~ 2 are 
over-estimated if derived from 24 yum flux densities, and discuss 
how this effect can be due to enhanced PAH em ission with resp ect 
to local templates (see also Elbaz et al. 2010; E lbaz et alj|201 ih . In 
NI2, the SED library from C harv & Elbaz qOolF was used, long- 
ward of 6 /im rest-frame, and the IRS spectra data were stacked to 
improve the faint signal of sources at z ~ 2. In this paper we adopt 
a different method to study the properties of z ~ 2 IR galaxies, as 
described in ^ a comparison of our results vs. those obtained by 
N12 is presented in 94.11 



Hereafter, we adopt the concordance cosmology {Hq = 
70 km sec ' Mpc ', n,„=0.3, and n^ =0.1. [Spergel et al.l I2OO3I) . 
Magnitudes are expressed in the Vega system. 



2 THE DATA 

The sample is composed of 24 luminous sources at z~2, selected 
by FIO in the Chandra Deep Field South (CDF-S), with 24 yum 
flux densities S(24 iim)~0.l4 - 0.55 mJy and at z = 1.75 - 2.40. 
The sample can be considered luminosity-selected, since all of the 
sources satisfying these selection criteria are considered. Given the 
redshift of the sources, their 24 pm flux densities translate roughly 
into infrared luminosity around lO'^Lo- We excluded two objects, 
originally defined as luminous IR galaxies (LIRGs; Ljj^ > 10" Lq) 
by FIO (L5511 and L62I1) and subsequently included in the 
ULIRG class thanks to the new IRS redshift measurements, since 
they would be the only two sources outside the HST ACS area 
(L551 1 is also outside the Herschel area). 

The sample benefits from a large amount of information avail- 
able for the CDF-S, from multi-band photometry to the recent 
Spitzer IRS spectroscopy. The IRS observations were performed 
at low-resolution in the observed 14-35 pm wavelength regime, i.e. 
sampling at z ~ 2 the important rest-frame PAH features at 6.2 pm 
and 7.7 pm, and the 9.7 pm silicate feature (see FIO). As reported 
above, in this work we also use recent far-IR data obtained from 
the Herschel satellite, which c onsists of PACS data at 70, 100 and 
160 pm from the PEP survey jLutz et alj|2oTTh . We use the PACS 
blind catalogue vl. 3 down to 3cr limits of 1.2, 1.2 and 2.0 mJy at 70, 
110 and 160 pm, respectively. Herschel data have been matched to 
24 pm Spitzer MIPS s ources jMagneUi et alj|2009l) u sing the like- 
lihoo d ratio technique ( Sutherland & Saunders! 19921; ICiliegi et alj 
l200lh . as described in iBerta et al.l fcoilh . Shorter wavelength in- 
formation has been included by matching the 24 yum sources to 
the multi-band (from UV to Spitzer IRAC bands) GOODS-MUSIC 
photometric catalog dSantini et al. I I2OO91V In particular, we cross- 
correlated the 24 pm selected ULIRGs with the PEP catalogue, 
considering the positions of the 24 pm sources already associated 
to the PEP ones. Among the 24 ULIRGs, 21 sources have counter- 
parts in PACS. For the 3 sources undetected by Herschel (U5050, 
U5152, and U5153), we consider the conservative 5cr upper lim- 
its (2.2, 2.0 and 3.0 mJy at 70, 100 and 160 pm, respectively). For 
PACS detections, the typical separation between 24 ytim and PEP 
sources is less than 1", with the exception of U5805 and U16526 
(~4"). By visually checking the far-IR images, these sources ap- 
pear as blends of more than one 24 ytim source. In these cases, the 
PACS flux densities are considered as upper limits. Photometry at 
16 pm, obtained by FIO from IRS data, has also been included. 

Finally, we consider the recently pu blished Chandra 4 Ms 
point-like source catalogue in the CDF-S jXue et alj|2oTll) . Eight 
of the 24 sources have an X-ray counterpart within 1 " . 

Table[T]lists the source names, redshifts from optical and IRS 
spectroscopy (see FIO), source IDs in the GOODS-MUSIC cata- 
logue, flux densities and associated errors from mid-IR (16 pm) to 
far-IR (160 pm). 



3 SED DECOMPOSITION 

The IR energy budget of a galaxy can be mainly ascribed to stellar 
photosphere and star-formation emissions and accretion processes; 
to estimate the relative importance of these three processes, a 
proper SED decomposition should be carried out. Disentangling the 
different contributions to the total SED is becoming more and more 
effective with the advent of the Spitzer and Herschel satellites. In 
this regard, many studies have been performed to compare the full 
range of observed photometric data with expectations from a host- 
galaxy component and models for the circum-nuclear dust emission 
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Table 1. The sample: multi-band information 



Name 


Zopt 


ZIRS 










S 100(1/11 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


U428 


-1.664 


1 783 


8053 


65 + 3 


2g9 + 5 




9070 + dsn 


< 3300 


U4367 


-1.762 


1 624 


3723 


101 + 8 


152 + 4 


< 1800 


< 1900 


7970 + S'^O 


U4451 


— 1 684 


1 875 


3689 




199 + 4 


< 1800 


7000 + ISO 


5110+ 550 


U4499 


-1.909 


1.956 


70361 


< 35 


167 + 3 


< 1800 


< 1900 


Mqn + SRO 


U4631 


1.896 


1 841 


7087 


< 21 


275 + 5 


< 1800 


< 1900 




U4639 


2.130 


2 112 


7553 


< 14 


213 + 4 


< 1800 


< 1900 


^570 + SSO 


U4642 


— 1 748 


1 898 


8217 


65 + 9 


242 + 5 


< 1800 


1 1 SO + ^ftO 


< 3300 


U4812 


1 910 


1 930 


13175 


< 53 


295 + 6 


9SQn + 410 


10000 + ^SO 


94R'^0 + 7'^0 


U4950 


2 291 


2 312 


15260 


453 + 8 


557 + 6 


< 1800 


70^iO + 490 


< 3300 


U4958 


2 145 


2 118 


15483 


J()9 + [Q 


232 + 4 


i^f;o+ 410 


7(^Q0 + ^SO 


^770 + SSO 




— 1 720 


1 938 


13758 


51+7 


197 + 4 


< 1800 


< 1900 


< 3300 


U5059 


-1.543 


1.769 


13887 


< 37 


272 ± 5 


1450 ± 380 


1750 ± 360 


4400 ± 610 


U5150 


-1.738 


1.898 


15083 


65 ± 11 


277 ±4 


< 1800 


4350 ± 440 


5690 ± 700 


U5152 


-1.888 


1.794 


70066 


66 ±9 


267 ± 8 


< 1800 


< 1900 


< 3300 


U5153 


-2.030 


2.442 


70054 


< 34 


166 ±3 


< 1800 


< 1900 


< 3300 


U5632 


1.998 


2.016 


9261 


83 ±2 


427 ± 5 


2290 ± 530 


4250 ± 360 


11470 ±560 


U5652 


1.616 


1.618 


6758 


154 ±9 


343 ±4 


1600 ± 380 


6920 ± 380 


17220 ± 740 


U5775 


-1.779 


1.897 


9361 


< 31 


164 ±4 


< 1800 


1440 ± 360 


3930 ± 670 


U5795 


-1.524 


1.703 


11201 


< 31 


250 ±4 


< 1800 


< 1900 


5030 ± 670 


U5801 


-1.642 


1.841 


12003 


< 35 


186 ±4 


< 1800 


< 1900 


2310 ±550 


U5805 


-2.093 


2.073 


12229 


68 ±6 


172 ±4 


< 1800 


1910''±360 


5100^ ±670 


U5829 


-1.597 


1.742 


9338 


57 ±5 


185 ±4 


< 1800 


1920 ± 380 


4320 ± 550 


U5877 


-1.708 


1.886 


13250 


179 ± 8 


364 ±5 


4600 ± 420 


7460 ± 440 


11570 ±870 


U16526 


-1.718 


1.749 


3420 


69 ±4 


306 ± 9 


4130° ±920 


8600'' ± 390 


15970" ±760 



Notes — (1) Source name (FIO); (2) optical redshift (negative and positive values are referred to photometric and spectro- 
scopic redshifts, respectively; see FIO for detai ls); (3) redshift derived from IRS spectroscopy (FIO); (4) ID from the v2 
GOODS-MUSIC c atalogue 1 Santini et alj|2009h : (5) 16 fim flux dens ity from FIO, in ^Jy; (6) 24 ^(m flux density in fi}y 
iMagneUi et al.l2009l) : (7), (8), (9) Herschel PEP flux density in //Jy <Lutz et al.ll201 ll) . " indicates a possible contribution 
from a nearby source. 



(e.g. JPolletta et tJj2008l ; lHatziminaoglou et alj2008l ; lvignali et alj 
I2OO9I ; IPozzi et al.ll201(]h . Here, together with the full-band photo- 
metric SED, we benefit also from the IRS spectra, t hat we combine 
with the photometric datapoints ( Sec. 13.21 see also Mullanev et alj 
l201lLlAlonso -Herrero et al.l20ll] for examples of combined photo- 
metric/spectroscopic data analysis). The IRS spectra provide an im- 
portant diagnostic, sampling the rest-frame mid-IR spectral range 
(~5 - 12 jim for our sample), where the difference between star- 
burst and AGN is strongest. In this wavelength range, starburst 
galaxies are generally characterized by prominent polycyclic aro- 
matic (PAH) features and weak 10 jim continuum, whereas AGN 
display weak or no PAH features plus a strong continuum (e.g. 
iLaurent et al]|2000l) . 



3.1 AGN and stellar components 

We have decomposed the observed SEDs using three distinct com- 
ponents: stars, having the bulk of the emission in the optical/near- 
IR; hot dust, mainly heated by UV/optical emission due to gas 
accreting onto the super-massive black hole (SMBH) and whose 
emission peaks between a few and a few tens of microns; cold dust , 
principally heated by star formation (w e refer to IPozzi et ^|201(]| : 
but see also lHatzimin aoglou et al.'2Q08' for a detailed description of 
the properties of the AGN and host-galaxy (stars -(-cold dust) com- 
ponents). Here we report only the most important issues concerning 
this analysis, with the AGN component being the main focus. 

The stellar component has been modelled as the sum of simple 
stellar population (SSP) models of solar metallicity and ages rang- 



ing from 1 Myr to 2.3 Gyr, which corresponds to the time elapsed 
between z = 6 (the redshift assumed for the initial star formation 
stars t o form) and z ~ 2 in the adopted cosmology. A Salpete^ 
( Il955h initial mass function (IMF), with mass in the range 0.15- 
120 Mq, is assumed. The SSP spectra have been weig hted by a 
Schmidt-like law of star formation (see lBerta et al.ll2004 



Tr - t Tn - t 

SFR{t) = — xexp - 



Tg 



Tar, 



(1) 



where Tq is the age of the galaxy (i.e. of the oldest SSP) and t,/ 
is the duration of the burst in units of the oldest SS P. A common 
value of extinction is applied to stars of all ages, and a lCalzetti et alJ 
|200(]|) attenuation law has been adopted (Rv-4.Q5). 

To account for emission above 24 yum, a component com- 
ing from colder, diffuse dust, likely heated by star-formation pro- 
cesses, has been included in the fitting procedure. It is represented 
by templates of well-studied starburst galaxies (i.e. Arp 220, M82, 
M83, NGC 1482, NGC 4102, NGC 5253 and NGC 7714) and 
five additi onal host-g alaxy average templates obtained recently by 
Mull anev et al] fcoili) from the starburst templates of IBrandl et al] 
(i2006l) . This set of five templates have been included since they 
properly reproduce the relative PAH strengths in the average IRS 
spectra of our ULIRG sample (see FIO). 

Regarding t he AGN compone nt, we have used the radiative 
transfer code of iFritz et al.l ( l2006h . This model follows the for- 
malism developed by diff e rent authors (e.g.. Pier & Krolik 19921 : 
iGranato & Danesd Il994t lEfstathiou & Rowan- Robinson Il995h . 
where the IR emission in AGN originates in dusty gas around 
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the SMBH with a "flared disk", "smooth distribution". Recently, 
this model has been widely used and found to successfully repro- 
duce the photometric data, including the 9.7 iim silicat e feature 
in em ission observed for type-I AGN (e.g., Siebenmo rgen et alj 
l2005h . Recent high-res olution, interfero metric mid-IR observations 
of nearby AGN (e.g., Ijaffe et al.|[200i) have confirmed the pres- 
ence of a geometrically thick, torus-like dust distribution on pc- 
scales; this torus is likely irregular or "clumpy". Over the last 
decade, many codes h ave been developed to deal with clumpy 
dust distributions (e.g..lNenkova et al.ll2002l : iNenkova et~^l2008l : 



iHonig & KishimotdllOim 



According to iDullemond & van Bemmell ( l2005h . the two 
models (smooth and clumpy) do not differ signi ficantly in re - 
producing sparse photometric datapoints (see also lElitzuJl2008h . 
The main difi^erence is in the strength of the silicate feature ob- 
served in absorption in objects seen edge-on, which is, on aver- 
age, weaker for clumpy models with the same global torus param- 
eters. The comparison between the two models has been applied 
only to few sampl es including IRS data and the results are not con- 
clusive. Recently, Ivignali et alj j201 ih made a comparison of the 
smooth vs. clumpy models for the matter responsible for repro- 
cessing the nuclear component of a hyper-luminous absorbed X-ray 
quasar at z~0.442 (IRAS 09104+4109), for which both photomet- 
ric and spectroscopic (IRS) data were available. While smooth so- 
lutions (the F06 model) are able to rep roduce the complete dataset, 
clumpy models jNenkova et alj[200s ) have problems in reproduc- 
in g the source photomet ry and spectroscopy at the same time. 
In iMullanev et 

aD j201lh . smooth vs. clumpy models are tested 
against a sample of ~10 local 5wi/f/BAT AGN with prominent 
emission at IRAS wavelengths. In this case, the authors claim that 
clumpy solutions reproduce the data better; smooth model parame- 
ters produce a much wider range of SED solutions, i.e. this model 
seems to produce too degenerate solutions. A further, more com- 
plete and extensive comparison of smooth vs. clumpy solutions is 
presented in Feltre et al. (submitted), where the theoretical SED 
shapes and the detailed spectral features of the two classes of mod- 
els (i.e. F06 for the "smooth distribution" and iNenkova et al ] |2008l 
for the "clumpy distribution") are compared using a large compi- 
lation of AGN with IRS spectroscopic data. Overall, similar results 
to those obtained by Dullemond & van Bemmel (2005) are derived, 
i.e. SED fitting applied to both photometric and spectroscopic data 
is not a sufficiently reliable tool to discriminate between the smooth 
and the clumpy distributions. We remind the reader that in the 
present paper we are focusing on the torus "global" energy out- 
put (i.e., the relevance of accretion-related emission with respect to 
the total source SED), not on the details of the torus structure and 
geometry, so the choice of the adopted model does not critically 
influence the results; as shown by Feltre et al. (submitted), the two 
models provide consistent results in terms of energetics. 

3.2 SED-fltting procedure 

In most previous work using the F06 code, only photometric dat- 
apoints were used, and the quality of the fitting solutions was 
estimated using a standard x~ minimization technique, where 
the observed values are the photometric flux densities (from 
optical-to-mid-IR/far-IR) and the model values are the "synthetic" 
flux densities obtained by convolving the sum of stars, AGN, 
and starburst components t hrou gh the filter respons e curves (see 
iHatziminaoglou et alj|200^ . In Vignali et alj ( l20Ilb . the spectro- 
scopic information was taken into account a posteriori in order to 
discriminate among different photometric best-fitting solutions. 



Here we propose a first attempt to simultaneously fit the pho- 
tometric and spectroscopic datapoints using smooth torus models 
by transforming the mid-IR 14-35 ;um observed-frame spectra into 
"narrow-band" photometric points of 1 yum band-width in the ob- 
served frame (i.e., subdividing the spectral transmission curve in 
sub-units), and estimating the coiTesponding fluxes and uncertain- 
ties using ordinary procedure of filter convolution and error prop- 
agation. The "new" filters have been chosen to achieve, in each 
wavelength bin, a sufficient signal-to-noise ratio without losing too 
much in terms of spectral resolution, which is needed to reproduce 
the 9.7 yum feature, when present. To take into account slit loss ef- 
fects, IRS data have been normalized to the 24 yum ffux density, 
deriving normalization factors between ~1.2 and 1.8. We would 
like to remark here that the F06 code does not necessarily impose 
the presence of an AGN component, i.e., solutions with only stellar 
emission are possible, as described in 34.11 Furthermore, the rel- 
ative normalization of the optical/near-IR component and the far- 
IR emission of the host galaxy is f ree, given the ex tremely com- 
plex physical relation of the two (i.e. lSaes et al.l2O10l) . Overall, the 
SED-fitting procedure ends with 1 1 free parameters: six are related 
to the AGN, two to the stellar component, and one to the starburst. 
The further two free parameters are the normalizations of the stel- 
lar and of the starburst components; the torus normalization is esti- 
mated by difference, i.e., it represents the scaling factor of the torus 
model capable to account for the data-to-model residuals once the 
stellar components have already been included. Here we brieffy re- 
call the parameters involved in our SED fitting analysis, and refer 
to F06 and Feltre et al. (submitted) for a detailed description of 
the AGN model parameter. The six parameters related to the torus 
are: the ratio R„,axlRmm between the outer and the inner radius of 
the torus (the latter being defined by the sublimation temperature 
of the dust grains); the torus opening angle 0; the optical depth r 
at 9.7 /im (t9 7^,„); the line of sight 6 with respect to the equato- 
rial plane; two further parameters, y and a, describe the law for the 
spatial distribution of the dust. In the currently adopted version (see 
Feltre et al., sub mitted), the "smoo th" torus database contains 2368 
AGN models. In lPozzi et alj j2010.) (see also lHatziminaoglou et alj 
( l2008h ). the degeneracies related to the six torus parameters are 
extensively described: in fact, various combinations of parameter 
values are equally able to provide good results in reproducing a set 
of observed data points. In particular, the optical depth T9 7^„, has 
the largest effect on the fit. The infrared luminosity, provided by 
the SED-fitting code, is robustly determined (within ~0.I dex) and 
appears "solid" against parameter degeneracies. 

Concerning the parameters associated with the other compo- 
nents, two are related to the stellar emission: t,/-, i.e. the parameter 
of the Schmidt-like law for the star formation, and the reddening 
E{B — V). Regarding the starburst component, the free parameter is 
related to the choice of the best-fitting template among the starburst 
library. 

The given number of free parameters means that the accept- 
able solutions, within l(3)cr confidence levels, are derived, for 
each source, by considering the parameter regions encompassing 
x\ji,j+(\2.65, 28 .5), respectively , in presence of II degrees of free- 
dom (d.o.f.; see lLampton et al.lll976) . 



4 AGN FRACTION 
4.1 SED-fltting results 

In Fig.Ia,b the observed UV-I60 /jm datapoints (filled red points) 
are reported along with the best-fitting solutions (black lines) and 
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Figure 1. a) Rest-frame broad-band datapoints (red dots) compared with the best-fit model obtained as the sum (solid black line) of a stellar (red dotted line), 
an AGN (blue dashed line) and a starburst component (green dot-dashed line). IRS spectra are shown as magenta lines. The area filled with diagonal lines 
represents the AGN solutions at the 3<t confidence level. 



the range of AGN models within the 3(T confidence level (filled 
region). All the sources need a host galaxy (red dotted-line) and 
a starburst component (green dot-dashed line). The host galaxy 
dominates the UV-8 jim photometry (at z ~ 2, the IRAC 8 /jm 
band samples the 2.7 fim rest-frame), while the starburst compo- 
nent dominates at longer wavelengths. For the three sources with 
no PACS detection (U5050, U5 152, and U5153), the starburst com- 
ponent is required by the SED-fitting procedure to reproduce the 
mid-IR spectral data, although its shape is not well constrained. 

Regarding the AGN component, our goal is to check whether 
its presence is required by the data and, in this case, to estimate its 



contribution to the IR luminosity. The SED-fitting procedure found 
that for all but one source (U428) the presence of an AGN is consis- 
tent with the photometry, although for only nine sources (~35% of 
the sample: U4639, U4950, U4958, U5150, U5152, U5153, U5652, 
U5805, and U5877) its presence is significant at the 3cr confidence 
level (i.e., solutions with no torus emission have x^^Xmi„+2S.5). 
However, in these sources as well, the AGN component far from 
dominates the whole spectral range, but emerges only in the narrow 
2-10 fjm range, where the stellar emission from the host galaxy has 
a minimum while the warm dust heated by the AGN manifests itself 
(e.g. . .Laurent et al...2000. : .Gruppioni et al...2010.) . 
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In Fig. |2] we report a zoom of the SED over the 2-20 yum 
range in order to visualize the AGN emission for all of the nine 
sources where such component is required at the 3cr confidence 
level. At a visual inspection, the presence of a nuclear component 
can be inferred for sources with a power-law SED (i.e., U4950 and 
U5877), for sources where the stellar component alone is not able 
to reproduce the datapoints around 3 ^m (i.e., U4958, U5153) or 
for sources where the starburst component, normalized to the far- 
IR datapoints, has already declined around 5-6 pm (i.e. U5152, 
U5153, and U5805). Finally, there are few cases (i.e., U4639) 
where visual inspection of the SED decomposition is not strongly 
suggesting an AGN component, which is however required by the 
SED-fitting analysis. 



The likely AGN in these sources is either obscured or of low- 
luminosity, although a combination of both effects is plausible as 
well. While our analysis is able to place constraints on the pres- 
ence of an AGN, we cannot draw any firm conclusion on either the 
obscured or the low-luminosity hypothesis for the AGN emission 
from mid-IR data. Despite the uncertainties affecting the estimate 
of th e gas column density A^^ derived from the dust optical depths 
(e.g., iMaiolino et J]|200lh . what we observe is that for all of the 
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Figure 2. Rest-frame SEDs and datapoints, as in Fig. [T] wliere a zoom in the 2-20 /jm wavelengtli range is shown for the nine sources where an AGN 
component is required at the 3 sigma confidence level from the SED-fitting analysis. 



sources with an AGN component, a certain level of obscuratiorQ 
(10^^ i,NH i 10^'' cm.-^) is required. 

In particular, the three sources where the presence of an AGN 
is particularly evident in the mid-IR regime from our SED de- 
composition (U4950, U4958, and U5877) - as already pointed out 
by FIO on the basis of the presence of a powerlaw mid-IR SED 
(U4950 and U5877), optical emission lines (U4958, with an appar- 
ently broad C ni] feature) and optically unresolved nucleus (U4950 
and U5877) - also show a relatively bright X-ray counterpart (see 
Sec. 14.21 ). The optical depth at 9.7 /jm of the best-fitting solutions 
corresponds to Nfj ~ (2 - 7)xl0"^ cm"- (i.e., in the Compton-thin 
regime; < N/j >= 6x10^^ cm"^ is obtained once the solutions at 
the 3o" level are considered). Of the remaining six sources, one 
(U4639) has an association with a relatively weak X-ray source, 
while for the others only an upper limit to the X-ray emission can 
be placed (see Table[2ll. For these six sources, the SED-fitting pro- 
cedure indicates an obscuration still in the Compton-thin regime. 



The column density has been derived from the the optical depth at 9.7/jm, 
using the Galactic e xtinction law (Cardelli et al. 1989) and dust-to-gas ratio 
iBohlinetal.1 19781) . 



but higher (up to 4x10^' cm"^, with < Nh >= 2x10^'' cm"^ when 
the solutions at the 3it level are taken into account) than for the 
previous three sources. 

Further insights on the properties of these sources, hence on 
the nature of their broad-band emission, will be discussed using 
X-ray diagnostics (see Sec. l4.2t . 

Turning now to source energetics, for the nine objects with 
an AGN component from the mid-IR, we computed the total and 
nuclear luminosity in three different spectral ranges: the whole (8- 
1000 jjm) IR range, the mid-IR (5-30 fim) range (partially sampled 
by IRS), and the narrow 2-6 yum wavelength interval. We find that 
the 8-1000 lum luminosity is always completely dominated by star- 
formation emission, the AGN nuclear contribution being i 5% 
(see Fig. |3] bottom panel), and that in only one source (U4950) 
out of the nine is the nuclear component contribution significant 
(~20%). Our finding that starburst processes dominate the 8-1000 
jim emission is consistent with the FIO and N12 conclusions. 

In the 5-30 yum range, where the re-processed emi ssion from 
the d ust surrounding the nuclear source peaks (e.g., ISilva et al] 
l2004h . we find a larger AGN contribution (i.e. ~25%; see Fig. [3] 
middle panel). As shown in Fig. la,b, at A~\Q pm the nuclear com- 
ponent typically starts being overwhelmed by the starburst emis- 
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z 



2.2 2.4 




ISOPHOT spectral observations, was ascribed to interstellar dust 
emission at temperatures of ~10'K). 

The dominance of the AGN emission in the narrow 2-6 yum in- 
terval is not in contrast with N12 conclusions, i.e. that these sources 
in the mid-IR are dominated by the PAH features, once the different 
spectral range is taken into account. 

In Table [3] we report the total IR luminosities (8-1000 yum), 
the AGN fractions over the three wavelengths discussed above (for 
the nine sources with an AGN emission detected at least at the 3cr 
confidence level), and both observed and predicted 2-10 keV lumi- 
nosities for our sample. The latter l uminosities h ave been derived 
from the 5.8 yum luminosity, using the lFiore et al.l (2009) correlation 
(see Sec. l4.2l for further details), and only for the nine sources with 
an AGN component detected at least at 3(T confidence (last col- 
umn of Table[3]l. The integration of the SED over the 8-1000 /jm 
range confirms the ULIRG classification {Lik > lO'^L©), inferred 
by FIO on the basis of the 24 yum flux densities and redshifts, 
for 14 out of the 24 sources, the remaining 10 sources showing 
slightly lower IR luminosities, between 0.5x10'' and 10'^ L©. The 
fact that 24 yum-based measurements tend to over-estimate the Lm 
for z~2 sources is consistent with other works based on Herschel- 
PEP data (e.g.,l Elbaz et al. 2010; N12) and stacking methods (e.g., 
IPapovich et alJuOOTl) . The final IR luminosity range of our sam- 
ple is 0.5-2.8x10'^ Lq (< l-m >= 1-4x10'- Lq, with a dispersion 
a- = 7x10" Lo). 




Figure 3. Fractional contribution due to the AGN component in the 2-6 jjm 
(top panel), 5-30 fim {middle panel), and 8-1000 /im (bottom panel) range. 
The error bars account for the AGN model dispersion at the 1 <t confidence 
level. Red points indicate the sources where an AGN component is detected 
at the 3o" confidence level from the SED fitting, and the triangles those with 
X-ray emission pointing clearly towards an AGN classification. 



sion. We note that the mid-IR AGN/starburst relative contribution 
is also discussed by FIO, adopting a completely different method 
than ours (i.e., scaling the SED of Mrk 231 to the 5.8 yum contin- 
uum and fitting th e residual SED with the average starburst from 
iBrandl et alj200 ^. Excluding from their analysis the three sources 
with the strongest AGN evidence from either X-ray or optical data 
(U4950, U4958, and U5877), FIO found a nuclear fraction of -20% 
(see their Fig. 22, top panel), which is consistent with our results. 

The only wavelength range where we find that the AGN over- 
comes the galaxy emission and contributes to ~60% of the emission 
is the narrow 2-6 yum interval (see Fig.|3] top panel). The power of 
the near-IR spectral range to detect obscured A GN emission con- 
, lRisalitietal]|2006l 2010 using L-band 



firms previous results (e.g. 
spectroscopy, but see also 
IR excess continuum emission, detected in disk galaxies thanks to 



4.2 X-ray results 

The 4 Ms X-ray source catalog in the CDF-S jXue et al.l201lh pro- 
vides additional information for the sub-sample of eight matched 
sources (see Table|2j. The depth of the Chandra mosaic in the field, 
though variable across its area, provides constraints down to very 
faint flux limits (~ 10"'^ erg cm"^ s"' in the 0.5-2 keV band). In 
the following, we use the X-ray information to provide an indepen- 
dent estimate of the AGN presence in our sample of IR-luminous 
galaxies and, for sources with most counts, characterize their X-ray 
emission. We also provide upper limits to the X-ray luininosity of 
the sources which are not detected in the 4 Ms CDF-S image^ 

A basic source classification is reported by Xue et al.! ( l20Ilb 
(AGN/galaxy/star), where five different classification criteria are 
adopted to separate AGN from galaxie^ at least one of these cri- 
teria should be satisfied. According to such classification, only two 
sources of the present sample are classified as galaxies, U428 and 
U4639, the remaining six being AGN. However, in the following 
we will provide indications for a "revision" of this classification 
using all of the available X-ray information (e.g., spectral proper- 
ties, luminosities, count distribution). 

In Table [2] the source classifications from X-ray data 
JXue et al]|201 ll and current work, respectively), along with the re- 
sults from the SED-fitting analysis (see Sec. l4.lt . are presented. 

X-ray luminosities, coupled with X-ray spectral analysis (see 
below), indicate the presence of an AGN in the three sources 
(U4950, U4958, and U5877) for which the SED fitting analysis 



Lu et where a weak near- 



^ The criteria to classify a source as an AGN are: high X-ray luminosity 
(above 3 x 10"*^ erg s"'); flat effective photon index (F < 1.0); X-ray-to- 
optical flux ratio log(fx/fR)> —1; an X-ray lumino sity at least three time s 
higher than that possibly due to star foiTnation (see [Alexander et al]|2005l) : 
presence of either broad or high-io nisation emission lines in the optical 
spectra. See §4.4 of lXue et al.l201 ll for details. 



The AGN content in luminous IR galaxies at z~2 9 



Table 2. X-ray properties of the sample of z ~ 2 IR-luminous galaxies 



Name 


XID Fo.5_8teV log(L2-10<rfv) log(L2-|0 tfV.fit) 


Nh 


Class. X (Xu et al. 11) 


Class. X (revised) Class. SED 


(1) 


(2) (3) (4) (5) 


(6) 


(7) 


(8) (9) 



Sources with an X-ray counterpart in the 4 Ms CDF-S catalog 



U428 


579 


0.07 


41.9 






Gal 


Gal 


Gal. 


U4639 


555 


0.07 


42.0 






Gal 


Gal 


AGN 


U4642 


437 


0.29 


42.5 


42.8 


4.5+*-^ X 10^2 


AGN 


AGN 


Gal. 


U4950 


351 


6.64 


43.9 


44.2 


1.2+0:3 X 1023 
1.9!'^^ X 102* 


AGN 


AGN 


AGN* 


U4958 


320 


0.31 


42.3 


43.7 


AGN 


AGN 


AGN* 


U5632 


552 


0.09 


42.1 




AGN 


Gal. 


Gal. 


U5775 


360 


0.05 


41.8 






AGN 


Gal. 


Gal. 


U5877 


278 


3.34 


43.0 


44.0 


6.Q^_l° X 1023 


AGN 


AGN 


AGN* 



Sources without X-ray detection in the 4 Ms catalog 



U4367 


< 


42. 


.3 


Gal. 


U4451 


< 


42. 


.2 


Gal. 


U4499 


< 


42. 


,1 


Gal. 


U4631 


< 


42. 


,0 


Gal. 


U4812 


< 


42. 


,4 


Gal. 


U5050 


< 


42. 


.2 


Gal. 


U5059 


< 


42. 


.0 


Gal. 


U5150 


< 


42. 


.4 


AGN 


U5152 


< 


42. 


,5 


AGN 


U5153 


< 


42. 


,6 


AGN 


U5652 


< 


41. 


,8 


AGN 


U5795 


< 


41. 


,8 


Gal. 


U5801 


< 


41. 


,9 


Gal. 


U5805 


< 


42. 


,1 


AGN 


U5829 


< 


41. 


,9 


Gal. 


U16526 


< 


42. 


.2 


Gal. 



Notes — (1) Source name (FIO); (2) XID from lXue etai]j201lh : (3) observed-frame 0.5-8 keV flux in units of 10"'^ erg cm"^ g-l jxue et alj2dTTI) : 
errors on the X-ray fluxes (hence on the luminosities) vary from less than 10% for the few sources with most counts up to ~ 40% for the X-ray faintest 
sources; (4) logarithm of the rest-frame 2-10 keV luminosity derived from either the 0.5-8 keV flux or the 0.5-8 keV sensitivity map (upper limits in 
the lower panel). The flux (or upper limit) is converted into a luminosity by assuming a powerlaw model with r=1.4; (5) logarithm of the rest-frame, 
absorption-corrected 2-10 keV luminosity, in units of erg s"'; (6) column density. Both (5) and (6) have been derived directly from X-ray spectral 
fitting, which has been limited to the four sources with the highest counting statistics (see Sect. l4.2lfor details). W e note that X-ray luminosities - 
columns (4) and (5) - are differe nt because of the different spectral modeling and assumptions (see lXue et al.l20l"l1 and Sect. 14.2) ; (7), (8), (9) source 
classifications from X-ray data l lXue et al.ll20l"ll and current work) and from SED-fitting, respectively. The term "AGN" in the source classification 
based on SED fitting (last column) indicates that the nuclear component is detected at the 3it confidence level. * indicates the presence of an AGN 
from optical data (U4958 - optical emission fines; U4950 and U5877 - optical morphology, see FIO). 



and FIO already suggested an AGN. Also U4642 has an X-ray lu- 
minosity typical of AGN emission. 

In comparison with the 2 Ms Chandra data jLuo et alj|200i) 
used by FIO, U4958 represents a new AGN detection, although the 
presence of an active nucleus was already inferred by the mid-IR 
featureless spectra and the optical spectrum, showing a broad C iii] 
line and strong N v and C iv emission lines. 

All of the X-ray matched sources have a full-band (0.5-8 keV) 
detection. One source (U4958) has an upper limit in the soft band 
(0.5-2 keV); this result is suggestive of heavy absorption, since 
Chandra has its highest effective area, and hence best sensitivity, in 
this energy range. Basic X-ray spectral analysis (using an absorbed 
powerla w with photon index fixed to 1.8, as typically observed 
in AGN: IPiconcelli et al.ll200^ . actually confirms the presence of 
strong, possibly Compton-thick obscuratior^ towards this source 



3 A source is called Compton-thick if its column density A'h > l/crj ~ 
1.5 X lO^"* cnT-, where crj is the Thomson cross section. 



(Nh = 1.9+^-^ X lO^'' cm-2). Four sources have hard-band (2-8 keV) 
upper limits': U428, U4639, U5632, and U5775, all of which being 
characterized by limited counting statistics (12-26 counts in the 
full band). Their X-ray luminosity (see Table I2FI is consistent with 
star-formation activity, although the presence of a low-luminosity 
AGN cannot be excluded. 

U4950 and U5877 are characterized by -1420 and -360 
counts in the 0.5-8 keV band, which allow a moderate-quality 
X-ray spectral analysis. Both sources can be fitted using a pow- 
erlaw model, but the flat photon index (-1 and -l-l, respectively) is 
indicative of absorption, which has been estimated to be 1.2+0 j ^ 
10^3 cnT- and 6.0+| x 10^3 cm"^^ respectively, after fixing the pho- 
ton index to 1.8. The same model applied to U4642 data provides 



For the four sources with 2-8 keV upper limits, the 2-10 keV luminosities 
are reported in Table|2]as these sources were detected in the hard band. We 
note, in fact, that the values reported in the table are derived from the 0.5- 
8 keV fluxes, where all of these sources are detected. 
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good results, although the derived column density is poorly con- 
strained (A^H = 4.5+5 J ^ 10^2 Apart from U4950, U4958, 
U5877 and U4642, for the remaining X-ray matched sources, the 
limited counting statistics prevent us from placing constraints on 
any possible obscuration. 

For the sources which were not detected in the 4 Ms CDF- 
S observations, we derived rest-frame 2-10 keV luminosity upper 
limits by converting the full-band sensitivity map (weighted over a 
region surrounding the source of interest of 4" radius to minimize 
possible spurious fluctuations in the map) using a powerlaw with 
r=1.4 (see Table|2]l- The assumption of r=1.8 would imply lumi- 
nosities higher by ~ 0. 1 - 0.2 dex. The derived X-ray upper limits 
place strong constraints on the further presence of luminous AGN 
in our sample, unless they are heavily obscured. 

Following FIG, for the nine sources with an AGN from the 
SED fitting, we estimated the rest-frame X-ray emission using the 
mid-IR emission of the AGN component as a proxy of the nuclear 
powe r (Table|3}. In particular, we use the L2_io*ci/-5.8 yum relation 
fere et aU2009, ; log Lz-iot.i/ = 43.57 + 0.72 x (l ogLs.g^m- 44.2), 
with luminosities expressed in erg s"'; see also iLutz et alj|2004h . 
which was calibrated using unobscured AGN in the CDF-S and 
COSMOS surveys. 

In Fig. |4] the predicted and measured X-ray luminosities are 
compared. For the three sources with clear AGN signatures (also 
from X-ray data), the luminosities have been corrected for the ab- 
sorption (see Tableland above), while for the remaining sources, 
no correction has been applied, due to the lack of constraints to any 
possible column density from the low-count X-ray spectra. Over- 
all, it appears that the three sources showing clear AGN signatures 
in the mid-IR band and in X-rays (U4950, U4958 and U5877) are 
broadly consistent with the 1:1 relation, considering the dispersion 
in the mid-IR/X-ray relation, the observational errors, and the un- 
certainties in deriving the intrinsic 2-10 keV luminosity properly. 
The correction accounting for the absorption seems to be correct 
also for the most heavily obscured of our X-ray sources, U4958. 

For the other six sources, the expected 2-10 keV luminosity 
is a factor of ~10-100 higher than the measured X-ray luminosity. 
Although caution is obviously needed in all the cases where mid-IR 
vs. X-ray luminosity coiTelations are used (see, e.g., the discussion 
inlyignali et al. 2010), the obscuration derived from the SED fit- 
ting for the remaining six sources (~ 2x10^^^ cm"-) and the ratio 
between the measured and predicted luminosities suggest the pres- 
ence of significant obscuration in these sources. Our conclusion is 
robust against the possible presence of absorption affecting also the 
Ls ifim values, since in this case the predicted X-ray luminosities 
should be considered as lower limits. Sensitive X-ray data over a 
larger bandpass would be needed to test the hypothesis of heavy 
obscuration in these sources. 



5 SUMMARY 

In this paper we have studied the multi-band properties of a sam- 
ple of IR luminous sources at z~2 in order to estimate the AGN 
contribution to their mid-IR and far-IR emission. The sample was 
selected by FIO at faint 24 flux densities (5(24yum)~0. 14 - 
0.55 mjy) and z = 1.75 - 2.40 to specifically target luminosities 
around 10'^ L©, i.e. sampling the knee of the IR luminosity func- 
tion. 

We have extended the previous analysis by taking advantage 
of new far-IR data recently obtained by the Herschel satellite as 
part of the guaranteed survey "PACS Evolutionary Probe" (PEP, 
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Figure 4. Predicted vs. measured 2-10 keV luminosities for the nine 
sources where an AGN component is inferred from the SED fitting. The tri- 
angles indicate the three sources with cleai' AGN X-ray emission; for these 
three sources only, the 2-10 keV luminosities have been corrected for ab- 
sorption derived from X-ray spectral analysis data (see Table|2)- The dashed 
diagonal lines indicate ratios of 1:1, 10:1 and 100:1 between the predicted 
and the measured X-ray luminosity (from right to left). 



iLutz et albOll ), and of the recently published 4 Ms Chandra data 
dXue et al.l201lh . The available photometry, coupled with IRS mid- 
IR data, have been used to reconstruct the broad-band SEDs of 
our IR-luminous galaxies. These SEDs have been modeled using 
a SED-fitting technique with three components, namely a stellar, 
an AGN, and a starburst component. The most up-to-date smooth 
torus model by F06 have been adopted for the AGN emission. The 
major results of the work can be summarized as follows. 

• SED fitting indicates that emission from the host galaxy in the 
optical/near-IR and star formation in the mid-IR/far-IR is required 
for all of the sources. Presence of an AGN component is consistent 
with the data for all but one source (U428), although only for 9 out 
of the 24 galaxies (~35% of the sample) is this emission significant 
at least at the 3cr level. 

• Of the sub-sample of nine sources that likely harbour an AGN 
according to the SED fitting, we find that their total (8-1000 fim) 
and mid-IR (5-30 fim) emissions are dominated by starburst pro- 
cesses, with the AGN-powered emission accounting for only ~ 5% 
and ~ 23% of the energy budget in these wavelength ranges, re- 
spectively. We find that the AGN radiation overcomes the stellar + 
starburst components only in the narrow 2-6 yum range, where it 
accounts for ~60% of the energy budget. In this wavelength range, 
stellar emission has significantly declined and emission from PAH 
features and starburst emission is not prominent yet. 

• For this same sub-sample, the gas column densities (derived 
by converting the dust optical depth at 9.7 yum obtained from 
the SED fitting) are indicative of a significant level of obscura- 
tion. In particular, three sources, also detected as relatively bright 
X-ray sources (U4950, U4958, and U5877, see below), have < 
Nm >= 6x10"^ cm"^ (considering all the solutions at the 3(T con- 
fidence level), while the remaining six sources have < A'^ >= 
2x10" cm--. 
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Table 3. IR and X-ray luminosities, and AGN fractions for the sample of z ~ 2 IR-luminous galaxies 
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Name 


ZIRS 




fAGN 
8-1000//m 


MGN 


fAGN 


log(L2-iofev) 




Class.SED 


(1) 


(2) 


(3) 


(4) 


(5) 


(6)"' 


(7) 


(8) 


(9) 


U428 


1.78 


12.06 


- 


- 


- 


41.9 


- 




U4367 


1.62 


11.70 


- 


- 


- 


- 


- 




U4451 


1.88 


12.00 


- 


- 


- 


- 


- 




U4499 


1.96 


12.25 


- 


- 


- 


- 


- 




U4631 


1.84 


11.89 


- 


- 


- 


- 


- 




U4639 


2.11 


12.08 


0.03 


0.16 


0.46 


42.0 


43.7 


AGN 


U4642 


1.90 


11.83 


- 


- 


- 


42.8 


- 




U4812 


1.93 


12.45 


- 


- 


- 


- 


- 




U4950 


2.31 


12.17 


0.22 


0.63 


0.91 


44.2 


44.4 


AGN 


U4958 


2.12 


12.10 


0.03 


0.11 


0.54 


43.7 


43.7 


AGN 


U5050 


1.73 


11.90 


- 


- 


- 


- 


- 




U5059 


1.77 


11.95 














U5150 


1.90 


12.24 


0.04 


0.26 


0.61 




43.8 


AGN 


U5152 


1.89 


11.94 


0.07 


0.26 


0.60 




43.7 


AGN 


U5153 


2.09 


12.26 


0.10 


0.51 


0.78 




43.8 


AGN 


U5632 


2.02 


12.39 








42.1 






U5652 


1.62 


12.37 


0.03 


0.20 


0.29 




43.7 


AGN 


U5775 


1.90 


11.98 








41.8 






U5795 


1.70 


11.93 














U5801 


1.84 


11.70 














U5805 


2.03 


12.12 


0.02 


0.13 


0.71 




43.7 


AGN 


U5829 


1.74 


11.90 














U5877 


1.89 


12.37 


0.03 


0.13 


0.68 


44.0 


43.9 


AGN 


U16526 


1.72 


12.38 














Notes — 


(1) Source name 


(FIO); (2) IRS redshift (FIO); (3) logarithm of the total (AGN+stellar) IR (8- 



1000 /jm) luminosity (in Lq); (4), (5), (6) AGN contribution to the 8-1000, 5-30 and 2-6 fjm luminosity, 
respectively, for the nine sources with an AGN component at the 3 cr level; (7) logarithm of the 2-10 keV 
luminosity (erg s"'; see Sec. l4.2l and Table|2); (8) logarithm of the 2-10 keV luminosity predicted from the 
L2- iokeV-Lj^um relat ion for sources with an AGN component detected at the 3cr confidence level (erg s"'; 
see iFiore et alj|2009l and Sec. 14. 2t : (9) Source classification from SED fitting. "AGN" indicates that the 
nuclear component is detected at the 3ct confidence level. All of the remaining sources are consistent with no 
significant AGN emission from the SED-fitting analysis. 



• X-ray analysis confiiin that three sources (U4950, U4958, and 
U5877) are actually powered by an AGN at short wavelengths, and 
that this AGN varies from being moderately (U4950 and U5877) to 
heavily obscured, possibly Compton thick (U4958). The X-ray lu- 
minosity of U4642 is also suggestive of moderately obscured AGN 
emission. The remaining four sources detected by Chandra have 
X-ray emission consistent with star-formation processes. 

• For the six sources where the AGN is required only by SED 
fitting (i.e., no strong AGN emission is observed in X-ray), we es- 
timate an intrinsic X-ray nuclear luminosity from the AGN contin- 
uum at 5.8 jim. The ratio (from 10 up to 100) between the predicted 
and the measured luminosities suggests the presence of significant 
obscuration in these sources. 
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